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ABSTRACT 

We present an analysis of high-resolution spectra of three core-helium-burning 'clump' 
stars and two first ascent giants in the open cluster IC 4651. Atmospheric parameters 
(T e ff, log g, vt, and [Fe/H]) were determined in our previous study by Carretta et al. 
(2004). In this study we present abundances of C, N, O and up to 24 other chemical 
elements. Abundances of carbon were derived using the C2 Swan (0,1) band head at 
5635.5 A. The wavelength interval 7980-8130 A, with strong CN features, was analysed 
in order to determine nitrogen abundances and 12 C/ 13 C isotope ratios. The oxygen 
abundances were determined from the [Ol] line at 6300 A. Compared with the Sun 
and other dwarf stars of the Galactic disk, mean abundances in the investigated clump 
stars suggest that carbon is depleted by about 0.3 dex, nitrogen is overabundant by 
about 0.2 dex and oxygen is close to solar. This has the effect of lowering the mean 
C/N ratio to 1.36 ± 0.11. The mean 12 C/ 13 C ratios are lowered to 16 ± 2. Other 
investigated chemical elements have abundance ratios close to the solar ones. 

Key words: stars: abundances - stars: atmospheres - stars: horizontal branch - open 
clusters and associations: individual: IC 4651. 
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1 INTRODUCTION 

This work is continuing our efforts in studying evolution- 
ary abundance alterations in evolved stars of open clusters 
(Tautvaisiene et al. 2000, 2005; Mikolaitis et al. 2010). Our 
main aim is to determine detailed elemental abundances of 
carbon, nitrogen and oxygen, and carbon isotope 12 C/ 13 C 
ratios in stars of open clusters in order to better understand 
reasons of abundance alterations caused by stellar evolution. 
Information on abundances of heavier chemical elements will 
be used for deriving the time evolution of abundances in 
the Galactic disk under the Bologna Open Cluster Chemical 
Evolution (BOCCE) study (Bragaglia & Tosi 2006, Carretta 
et al. 2007, and references therein). In this work, our target 
of investigations is the open cluster IC 4651. 

The open cluster IC 4651 is an intermediate-age 
(1.7 Gyr) open cluster located 140 pc below the Galactic 
plane and 7.1 kpc form the Galactic centre (Meibom et al. 
2002; Pasquini et al. 2004). Meibom et al. (2002) provided 
calculations of the space motion and the Galactic orbit of the 



* Based on observations collected at ESO telescopes under pro- 
grammes 65.N-0286 and in part 169.D-0473 
f E-mail: sarunas.mikolaitis@tfai.vu.lt 



cluster. The orbital eccentricity is e = 0.19 and the mean ra- 
dius of galactocentric orbit is 8.6 kpc, its maximum distance 
from the Galactic plane is 190 pc (0:2000 = 17 ,l 24.8 m , foooo = 
-49°56.0';Z = 340.088°, b= -07.907°). 

Results of extensive photometric studies were published 
first by Eggen (1971) and Lindoff (1972) and later on by 
Anthony-Twarog & Twarog (1987, 2000), Anthony-Tworag 
et al. (1988). However, the most recent photometric study 
was performed by Meibom (2000) and Meibom et al. (2002). 
They combined photometric observations and radial veloc- 
ity measurements for the 44 single member stars down to 
V = 14.5 mag and determined E(B — V) = 0.10 mag, the 
distance d — 1.01 ± 0.05 kpc, and the mean radial velocity 
equal to -30.76 ± 0.20 kms" 1 . It was found that 37% of 
giant members are spectroscopic binaries with periods up 
to 5000 days, and 52% of the main-sequence and turn-off 
members are binaries with periods less than 1000 days. The 
estimated total mass of IC 4651 is « 630Mq (Meibom et al. 
2002). The turn-off mass of the IC4651 stars M = 1.69M 
was obtained by Carretta et al. (2004) reading the turn-off 
values on the Girardi et al. (2000) isochrones for solar metal- 
licity at the age of the cluster of 1.7 Gyr as determined by 
Meibom et al. (2002). 
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Figure 1. The colour-magnitude diagram of the open cluster 
IC4651. The stars investigated in this work are indicated by 
the filled squares. The stars of two other high-resolution spec- 
tral abundance studies are shown in this plot as well: the work 
of Pace et al. (2008) is marked by triangles and of Pasquini et 
al. (2004) - by diamonds. The diagram is based on Stromgren 
photometry by Anthony-Twarog & Twarog (2000). 



There were several photometric studies that evaluated 
the metallicity of IC4651. Based on uvby — Hp photome- 
try, [Fe/H] = 0.23 ± 0.02 was found by Anthony-Twarog 
& Twarog (1987), [Fe/H] = 0.18 ± 0.05 by Nissen (1988), 
and [Fe/H] = 0.077 ± 0.012 by Anthony-Twarog & Twarog 
(2000). 

High resolution spectroscopic data started to appear in 
the beginning of the millennium. Bragaglia et al. (2001) de- 
termined the mean metallicity of five evolved stars [Fe/H] = 
0.16 ± 0.01. Pasquini et al. (2004) provided the cluster 
metallicity [Fe/H] = 0.10 ± 0.03 from the analysis of 22 
faint main sequence stars of the cluster. Abundances of 
the iron peak, a-elements and lithium were investigated in 
their study as well. Carretta et al. (2004) found the average 
[Fe/H] = 0.11 ±0.01 for 5 evolved stars of the cluster, which 
we analyse further in this our work. Pace at al. (2008) pro- 
vided the abundance measurements of Fe, Ca, Na, Ni, Ti, 
Al, Cr, Si for 20 solar-type stars belonging to IC 4651 and 
found [Fe/H] = 0.12 ±0.05. And finally, Santos et al. (2009) 
derived the mean metallicity [Fe/H] = 0.15T0.02 for IC4651 
from five dwarf stars. 

In our work for the open cluster IC4651, the detailed 
abundance analysis of almost 30 chemical elements is done. 
Abundances of such key chemical elements as 12 C, 13 C, N, 
O as well as representatives of s- and r-processes are deter- 
mined. 

The colour-magnitude diagramme of IC4651 with the 
stars analysed in three most comprehensive chemical abun- 
dance studies (Pasquini et al. 2004; Pace et al. 2008 and this 
work) indicated is presented in Fig. 1. 




Figure 2. Field of 8 x 8 arcmin 2 centered on IC4651, with the 
programme stars indicated by their numbers according to Lindoff 
(1972). 



2 OBSERVATIONS AND METHOD OF 
ANALYSIS 

The spectra of five cluster stars (IC4651 27, 56, 72, 76 and 
146) were obtained with the spectrograph FEROS (Fiber- 
fed Extended Range Optical Spectrograph) mounted at the 
1.5 m telescope in La Silla (Chile). The resolving power is 
R = 48 000 and the wavelength range is AA 3700-8600 A. 
Three stars (27, 76, 146) belong to the red clump of the 
cluster, the IC 4651 72 star is a first-ascent giant, and the 
star 56 is an RGB-tip giant (see Fig. 1). The finding chart 
of the investigated stars is shown in Fig. 2. The log of ob- 
servations and S/N are presented in the paper by Carretta 
at al. (2004). 

In the same paper by Carretta et al. (2004), all the main 
atmospheric parameters for the observed stars were deter- 
mined. For the convenience we present them in this paper 
as well (Table 1). The effective temperatures were derived 
by minimizing the slope of the abundances from neutral Fe I 
lines with respect to the excitation potential. Using the line- 
depth ratios (LDR) technique, Biazzo et al. (2007) have 
determined higher effective temperatures for other clump 
stars in NGC4651 and have raised doubts that the temper- 
ature determinations by Carretta et al. were likely too low. 
We have checked dependences of the other chemical element 
lines with respect to the excitation potential and did not find 
slopes. Thus, we do not doubt in the correctness of effective 
temperature determinations for IC4651 giants by Carretta 
et al. 

The gravities (log g) were derived by Carretta et al. 
from the iron ionization equilibrium. In our study we found 
a very good agreement between neutral and singly ionised 
species of Cr and Ti, which strongly support the reliabil- 
ity of the atmospheric parameters, in particular the gravity 
values derived from the ionization equilibrium of Fe. The mi- 
croturbulent velocities were determined assuming a relation 
between log g and Vt ■ The ATLAS models with overshooting 
were used for the analysis. The Fel lines were restricted to 
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Table 3. Abundances relative to hydrogen [El/H]. The quoted errors, a, are the standard deviations in the mean value due to the line-to-line scatter within the species. The number of 
lines used is indicated by n. The last two columns give the mean [El/Fc] and standard deviations for the cluster stars 27, 72, 76 and 146. 



Species 




27 
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72 
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Mean 
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(j 


Tl 


fEl/Hl 


(j 


TL 


[El/Fc] 




C(C 2 ) 


-0.15 




1 


-0.64 




1 


-0.12 




2 


-0.17 




1 


-0.15 




1 


-0.26 


0.02 


N (CN) 


0.36 


0.07 


24 


-0.09 


0.11 


18 


0.35 


0.07 


24 


0.30 


0.06 


24 


0.27 


0.09 


23 


0.21 


0.04 


0([0 1]) 


0.05 
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0.09 
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0.01 


0.06 
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0.09 


0.05 
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-0.25 


0.02 
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0.04 


3 


0.10 


0.05 
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0.07 


0.04 
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0.01 
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0.05 


0.07 
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0.08 
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0.05 


0.09 
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0.05 
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0.03 
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0.05 
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0.07 


4 


0.13 
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0.25 


0.08 
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0.17 


0.02 
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0.05 


9 


0.14 


0.06 
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0.15 


0.06 


9 
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0.05 
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0.05 


0.03 


Ti I 


0.25 


0.08 


24 
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0.07 


9 


0.23 


0.09 


24 


0.18 


0.08 


26 


0.21 


0.10 


28 


0.11 


0.03 


Tin 


0.21 


0.09 
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0.09 
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0.07 
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0.09 


11 
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0.09 


11 


0.08 


0.04 
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0.20 


0.07 
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0.05 
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0.10 
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19 
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0.09 


24 
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24 
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0.05 


0.04 
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0.01 


6 


0.21 
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2.1 Estimation of uncertainties 

The sources of uncertainty were described in detail in our 
Paper I. 

The sensitivity of the abundance estimates to changes in 
the atmospheric parameters by the assumed errors (± 100 K 
for T e g, ±0.3 dex for log g and ±0.3 km s _1 for v t ) is illus- 
trated for the star IC4651 72 (Table 2). It is seen that possi- 
ble parameter errors do not affect the abundances seriously; 
the element-to-iron ratios, which we use in our discussion, 
are even less sensitive. The sensitivity of iron abundances to 
stellar atmospheric parameters were described in Carretta 
et al. (2004). 

The scatter of the deduced line abundances a, presented 
in Table 3, gives an estimate of the uncertainty due to the 
random errors, e.g. in continuum placement and the line pa- 
rameters (the mean value of a is 0.06). Thus the uncertain- 
ties in the derived abundances that are the result of random 
errors amount to approximately this value. 

Since abundances of C, N and O are bound together by 
the molecular equilibrium in the stellar atmosphere, we have 
also investigated how an error in one of them typically affects 
the abundance determination of another. A[0/H] = 0.10 
causes A[C/H] = 0.05 and A[N/H] = -0.10, A[C/H] = 0.10 
causes A[N/H] = -0.15 and A[0/H] = 0.02, A[N/H] = 0.10 
has no effect on either the carbon or the oxygen abundances. 



3 RESULTS AND DISCUSSION 

The abundances relative to hydrogen [E1/H0 and a (the 
line-to-line scatter) derived for up to 27 neutral and ionized 
species (including 13 C) for the programme stars are listed 
in Table 3. The average cluster abundances [El/Fe] and dis- 
persions about the mean values for IC4651 are presented in 
Table 3 as well. They are calculated from the results deter- 
mined for the stars 27, 72, 76 and 146. Due to the differ- 
ent [Fc/H], which is by 0.4 dex lower than of other cluster 
stars, the star 56 was not used in the average calculations 
even though its values do not change the average abundances 
much. For the majority of the chemical elements the changes 
are just ±0.01 — 0.02 dex, only for Si I and Sell the differ- 
ence is 0.04 dex, and for Vl - 0.07 dex. From its element 
to iron ratios, C/N and 12 C/ 13 C ratios, this star within er- 
rors of uncertainties is indistinguishable from other evolved 
stars of this cluster. However, we think that with more ade- 
quate model atmospheres this star can be investigated more 
accurately. In Fig. 4 we show an example of spectrum syn- 
theses for the [Ol] line in IC4651 56. In Fig. 3 a sample of 
its spectrum is shown as well. 

In IC4651, the majority of investigated chemical ele- 
ments have abundance ratios close to the solar ones. The 
mean cluster [a/Fe] = ±([Mg/Fe] ± [Si/Fe] ± [Ca/Fe] ± 
[Ti/Fe]) = 0.06 ± 0.07 (s.d.), which is close to the solar 
value. 

Abundances of Nal, All, Si I, Cal, Til and Nil were 
determined for the main-sequence stars of IC 4651 by Pace 
et al. (2008). The mean [El/Fe] ratios in these stars are very 
close to solar as well. 

1 In this paper we use the customary spectroscopic notation 
[X/Y]=log 10 (Ar x /iV Y ) 

star — l°glo(-^x/-^Y)© 
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Figure 4. Fit to the forbidden [O i] line at 6300 A in IC 4651 56. 
The observed spectrum is shown as a solid line with black dots. 
Synthetic spectra with [0/Fe]= 0.3, 0.08, and 0.13 are shown as 
solid gray lines. 



Table 4. The mean [El/Fe] for giant (G) and main-sequence (MS) 
stars of IC4651 investigated in this work, Pasquini et al. (2004) 
and Pace et al. (2008). 





This work 


Pasquini 


Pace 


Species 


[El/Fe] 


[El/Fe] 


[El/Fe] 


[El/Fe] 




G 


G 


MS 


MS 


Nal 


0.00 


0.19 


-0.09 


-0.03 


Mgi 


-0.05 


0.09 


0.13 




All 


0.03 


0.07 


-0.07 


-0.10 


Si i 


0.09 


0.08 


0.07 


-0.02 


Cai 


0.05 


0.00 


0.04 


0.04 


Sen 


0.01 


0.11 


-0.11 




Til 


0.11 


0.12 


0.08 


-0.02 


Tin 


0.07 


0.18 


0.00 




Cri 


-0.05 


-0.02 


0.11 




Nil 


0.07 


0.10 


0.01 


-0.02 



Pasquini et al. (2004) investigated spectra of both giants 
and main-sequence stars in IC4651. The authors expressed 
their strong believe that [Na/Fe] ratio is comprehensively 
higher in the giants in comparison to the main-sequence 
stars and that this is due to internal nucleosynthesis and 
mixing. However, neither Pace et al. (2008), neither our 
study may confirm this statement. In Table 3, we present the 
mean [El/Fe] for giant and main-sequence stars of IC4651 
investigated in our work, Pasquini et al. (2004) and Pace et 
al. (2008). In our work, the abundances of Na and Mg were 
determined with NLTE taken into account, and we do not 
find an overabundance of these chemical elements. 

In Fig. 4 of Paper I, we presented the radial distribution 
of some elemental abundances for BOCCE clusters analysed 
so far, and for others in recent studies. IC 4651 agrees well 
with results of other open clusters at the same R gc of 7.1 kpc. 

In the following sections we will discuss in more detail 
the results of carbon and nitrogen abundance determina- 
tions. 
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Figure 5. Small region of IC4651 27 spectrum (solid black line 
with black dots) at C2 Swan (0,1) band head 5635.5 A, plotted 
together with synthetic spectra with [C/Fe] values lowered by 
—0.2 dex (lower gray line), —0.25 dex (middle gray line) and —0.3 
(upper gray line). 



The mean nitrogen to iron abundance ratio in IC 4651 is 
[N/Fe] = 0.21 ± 0.04. This shows that nitrogen is overabun- 
dant in these evolved stars of IC4651, since [N/Fe] values 
in the Galactic main-sequence stars are about solar at the 
solar metallicity (c.f. Shi et al. 2002). Unfortunately, nei- 
ther carbon, neither nitrogen abundances were investigated 
in the main sequence stars of IC 4651 by Pace et al. (2008) 
and Pasquini et al. (2004). 

The mean C/N ratios in IC4651 is equal to 1.36 ±0.11. 
The smallest value of C/N = 1.12 was obtained for the star 
IC 4651 56. 

The 12 C/ 13 C ratios were determined for all programme 
stars from the (2,0) 13 C 14 N feature at 8004.728 A. In Fig. 6 
we show a small region of IC 4651 72 spectrum together with 
spectral syntheses obtained with three different values of the 
carbon isotopic ratio. We find that the mean 12 C/ 13 C ratios 
are about 16 ± 2 in the evolved stars investigated. 

The solar carbon and nitrogen abundances used in our 
work arc log^c = 8.52 and log^N = 7.92 (Grevesse & 
Sauval 2000), so the solar C/N = 3.98. The 12 C/ 13 C ra- 
tio in the solar photosphere is equal to 89 (Coplen et al. 
2002). 



0.75 



- 1 1 1 r 




8003.0 



8003.5 



8004.0 

X{k) 



8004.5 



8005.0 



Figure 6. Small region of IC4651 72 spectrum (solid black line 
with black dots) with 13 C 14 N feature. Grey lines show synthetic 
spectra with 12 C/ 13 C ratios equal to 13 (lower line), 15 (middle 
line) and 17(upper line). 



3.1 Carbon and nitrogen abundances 

The average value of carbon to iron ratio in IC 4651 is 
[C/Fe] = -0.27 ± 0.02. In Fig. 5, a fit to the IC4651 27 
spectrum at C2 5635.5 A is shown. 

We compared the carbon abundance in IC 4651 with 
carbon abundances determined for dwarf stars in the Galac- 
tic disk. Shi et al. (2002) performed an abundance analysis of 
carbon for a sample of 90 F and G type main-sequence disk 
stars using Cl and [Cl] lines and found [C/Fe] to be about 
solar at the solar metallicity. Roughly solar carbon abun- 
dances were found by Gustafsson et al. (1999) who analysed 
a sample of 80 late F and early G type dwarfs using the for- 
bidden [Cl] line. The ratios of [C/Fc] in our stars lie about 
0.3 dex below the values obtained for dwarf stars of the 
Galactic disk. 



3.2 Comparison of 12 C/ 13 C and C/N ratios with 
theoretical models 

The carbon and nitrogen abundances, C/N and especially 
the carbon isotope ratios 12 C/ 13 C are key tools for stel- 
lar evolution studies. Investigations of abundances of these 
chemical elements in atmospheres of clump stars of open 
clusters may provide a comprehensive information on chem- 
ical composition changes. The clump stars have accumulated 
all chemical composition changes that have happened dur- 
ing their evolution along the giant branch and during the 
helium flash. 

In Fig. 7 and 8, we compare the mean carbon isotope 
and C/N ratios of clump stars in different open clusters as 
a function of turn-off mass with the theoretical models of 
the 1 st dredge-up, thermohaline mixing (TH), thermohaline 
mixing together with rotation-induced mixing for stars at 
the zero age main sequence (ZAMS) having rotational ve- 
locities of 110 kms -1 , 250 kms -1 and 300 kms -1 computed 
by Charbonnel & Lagarde (2010); and Cool Bottom Pro- 
cessing model (CBP) by Boothroyd & Sackman (1999). 

The most resent modelling of extra-mixing processes 
was done by Charbonnel & Lagarde (2010). They are based 
on ideas of Eggleton et al. (2006) and Charbonnel & Zahn 
(2007) . Eggleton et al. (2006) found a mean molecular weight 
(fi) inversion in their IMq stellar evolution model, occurring 
after the so-called luminosity bump on the red giant branch, 
when the H-burning shell source enters the chemically homo- 
geneous part of the envelope. The /i-inversion is produced by 
the reaction 3 He( 3 He, 2p) 4 He, as predicted by Ulrich (1972). 
It does not occur earlier, because the magnitude of the p- 
inversion is small and negligible compared to a stabilizing 
/x-stratification. Following Eggleton et al. (2006), Charbon- 
nel & Zahn (2007) have computed stellar models including 
the prescription by Ulrich (1972) and extended them to the 
case of a non-perfect gas for the turbulent diffusivity pro- 
duced by that instability in a stellar radiative zone. They 
found that a double diffusive instability referred to as ther- 
mohaline convection, which has been discussed long ago in 
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6 - Mikolaitis et al. (2010); 7 - This work. 
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the literature (Stern 1960), is important in evolution of red 
giants. This mixing connects the convective envelope with 
the external wing of hydrogen burning shell and induces sur- 
face abundance modifications in red giant stars. 

Charbonnel & Lagarde (2010) also computed the mod- 
els of rotation-induced mixing for stars at the zero age 
main sequence having rotational velocities of llOkms -1 , 
250kms _1 and 300 km s" 1 . Typical initial ZAMS rotation 
velocities were chosen depending on the stellar mass based 
on observed rotation distributions in young open clusters 
(Gaige 1993). The convective envelpoe was supposed to ro- 
tate as a solid body through the evolution. The transport 
coefficients for chemicals associated to thermohaline and 
rotation-induced mixings were simply added in the diffu- 
sion equation and the possible interactions between the two 
mechanisms were not considered. The rotation-induced mix- 
ing modifies the internal chemical structure of main sequence 
stars, although its signatures are revealed only later in the 
stellar evolution. 

The models by Boothroyd & Sackmann (1999) include 
the deep circulation mixing below the base of the standard 
convective envelope, and the consequent "cool bottom pro- 
cessing" (CBP) of CNO isotopes. 

The theoretical models were compared to the observa- 
tional data of 12 C/ 13 C and C/N listed in Table 5, which 
we collected for clump stars of open clusters investigated by 
Mikolaitis et al. (2010), Smiljanic et al. (2009), Tautvaisiene 
et al. (2000, 2005), Luck (1994) and Gilroy (1989). From 
Gilroy (1989) we selected 4 clusters with well defined red 
clump stars. Luck (1994) derived carbon isotope ratios for 8 
open clusters, however only one cluster was included to our 
comparison since for other clusters it was very difficult to 
identify stars of red clump. The turn-off masses, ages and 
galactocentric distances were choosen from the most recent 
studies and used for displaying of other 12 C/ 13 C and C/N 
investigations for the same cluster, if available. 

In Fig. 7 and 8 we can see that for clusters with stars of 
smaller turn-off masses the 12 C/ 13 C and C/N values arc m a 
good agreement with both models of extra-mixing used for 
the comparison. However, 12 C/ 13 C values in the clump stars 
with turn-off masses of 2-3 Mq in most of the investigated 
clusters are lower than predicted by the available models 
and need modelling of larger extra-mixing. 

3.3 Final remarks 

Carbon and nitrogen are important products of nucleosyn- 
thesis processes in stellar interiors, and the evidence of their 
abundance variation during stellar evolution is a signature 
of physical mixing processes between the atmosphere and 
deeper layers of a star. Such abundance alterations may be 
well traced in open clusters. They provide a unique possibil- 
ity for investigation of a number of stars of nearly the same 
age, distance and origin, as open cluster stars are claimed 
to be formed in the same protocloud of gas and dust. Open 
clusters have a high reliability of mass, distance, evolution- 
ary phase and abundance determinations. 

Extra-mixing processes may become efficient on the red 
giant branch when stars reach the so-called RGB bump, and 
may modify the surface abundances. It is known that alter- 
ations of 12 C/ 13 C and 12 C/ 14 N ratios depend on stellar evo- 
lutionary stage, mass and metallicity (see Charbonnel et al. 
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Figure 7. The average carbon isotope ratios in clump stars of 
open clusters as a function of stellar turn-off mass. The result of 
this work is marked by the filled square; from Mikolaitis et al. 
(2010) and Tautvaisiene et al. (2000, 2005) - open squares; from 
Smiljanic et al. (2009) - open triangles; from Luck (1994) - re- 
versed open triangle; from Gilroy (1989) - open circles. The mod- 
els of the 1 st dredge-up, thermohaline mixing (TH) and rotation- 
induced mixing (V) arc taken from Charbonnel & Lagarde (2010). 
The CBP model of extra-mixing is taken from Boothroyd & Sack- 
mann (1999). 
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Figure 8. The average carbon to nitrogen ratios in clump stars of 
open clusters as a function of stellar turn-off mass. The meaning 
of symbols are as in Fig. 7. 



1998, Gratton et al. 2000, Chaname et al. 2005, Cantiello & 
Langer 2010; Charbonnel & Lagarde 2010 for more discus- 
sion). 

The comparison of the observational data with theo- 
retical models of stellar evolution shows that processes of 
extra-mixing in stars of open clusters with turn-off masses 
of 2-3 Mq are larger than predicted. 
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